Note and calculations concerning elastic dilatancy in 2D glass-glass liquid foams 
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When deformed, liquid foams tend to raise their liquid contents like immersed granular materials, 
a phenomenon called dilatancy. We have aready described a geometrical interpretation of elastic 
dilatancy in 3D foams and in very dry foams squeezed between two solid plates (2D GG foams). 
Here, we complement this work in the regime of less dry 2D GG foams. In particular, we highlight 
the relatively strong dilatancy effects expected in the regime where we have predicted rapid Plateau 
border variations. 

PACS numbers: 47.20.Dr, 83.80.Iz, 47.57.Bc, 68.03.Cd 



I. DILATANCY IN FOAMS 

Liquid foams (T), Q and granular materials both exhibit 
"dilatancy", described by Reynolds Q in the context of 
granular materials: upon deformation, because grains are 
forced to move while avoiding each other, the medium 
swells to some extent. In other words, the fluid volume 
fraction </> is increased (see Fig. [IJ . This effect can remain 
unnoticed in air. By contrast, a spectacular absorption 
of liquid Q is obtained upon deformation of an immersed 
granular sample. 

Because bubbles can deform individually, a foam might 




FIG. 1: Dilatancy in granular media and in foams. When 
a granular material (left), initially in a compact state, is de- 
formed, the steric interactions between grains cause the liq- 
uid content to rise. By contrast, in a foam (right), because 
the bubbles are deformable, the liquid volume fraction may 
remain constant. Static dilatancy means that if the liquid 
fraction tends to increase due to the applied deformation. 
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deform substantially without altering its liquid fraction 
(see Fig.QJ. Hence, there is a priori no reason why a foam 
should display dilatancy. In fact, dilatancy does exist in 
foams. This has been shown both experimentally [f| in 
the case of "dynamic" dilatancy (caused by a continuous 
foam deformation) and theoretically [6 8] in the case of 
"clastic" dilatancy (caused by a constant deformation of 
the foam). 

Dilatancy can be defined @ in terms of the osmotic 
pressure of a foam confined in a container, which corre- 
sponds to the force that must be applied externally to 
one of the confining walls if the latter is permeable to 
the liquid but not to the bubbles. The static dilatancy 
coefficient \ reflects the fact that the foam osmotic pres- 
sure varies with the deformation e when the fluid volume 
fraction <p is kept constant: 



Surprisingly, this coefficient can be either positive like 
in granular materials or negative [(|. In the case of a 
negative dilatancy coefficient, deforming the foam results 
in a tendency to expcll liquid and make the foam dryer. 

We have shown that the origin of this change in sign 
can be traced back to two different physical contribu- 
tions Q. On the one hand, deforming the foam implies 
an increase in the total of all Plateau border lengths in 
the sample (in the case of a 2D GG foam Q, these are 
the pseudo Plateau borders, i.e., those along the solid 
plates). As a result of this increase in length, when the 
total amount of liquid is kept constant, the Plateau bor- 
ders shrink, which raises the pressure difference between 
the gas and the liquid, resulting in an increase of the 
osmotic pressure, hence a positive contribution to dila- 
tancy. On the other hand, deforming the foam also im- 
plies an increase in the total surface area of the gas-liquid 
interfaces. Because the interfaces contribute negatively 
to the stress tensor in the foam (tensile contribution). 
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this has a negative contribution to the osmotic pressure, 
and hence, to dilatancy. Because the increase in total 
interfacial energy is directly related to the elasticity of 
the foam, this negative contribution to dilatancy coeffi- 
cient is proportional to the elastic modulus, as shown by 
Weaire and Hutzler Jf|, 

In the present work, with this geometrical interpreta- 
tion in mind, we conduct the explicit calculation of the 
clastic dilatancy of not too dry 2D GG foams (regimes 
A, B, C and D). In particular, we highlight the relatively 
strong dilatancy effects expected in the regime where we 
have predicted rapid Plateau border variations. 




II. GEOMETRY OF 2D GG FOAMS: FLOOR 
TILE VERSUS PANCAKE REGIME 

Let us first choose notations || [!(} to describe the ge- 
ometry of two-dimensional foams squeezed between two 
glass plates, called "GG" foams. 

We call L the perimeter of the bubble, defined as the 
perimeter of the rounded polygon that constitutes the 
bubble, as seen from above, and R the corresponding 
radius of curvature of the Plateau border (see Fig. [3J 
left). We call H the distance between both solid plates, 
and R ps the radius of curvature of the pseudo Plateau 
borders (see Fig. [5J bottom right). We also call the 
average bubble volume, and f2ii q the average volume of 
liquid per bubble. The liquid volume fraction <fi thus 
verifies: 



liq 



filiq = Q 



1 - 



(2) 



We assume that f2 remains constant: the applied stresses 
are not sufficient to compress the gas phase significantly 
unless the bubble size is on the order of a micron. 

The main contributions to the quantity of liquid per 
bubble are pictured on Fig. [2] and can be calculated from 
simple geometrical arguments: 



fi liq ~(2-7r/2)L^ s 



(2\/3 - tt) R 2 H (3) 



The first term corresponds to the pseudo Plateau bor- 
ders, which make the junction between the intcrbubblc 
films and the solid plates. As seen from above, they cor- 
respond to the light grey regions in Fig. [2j Each portion 
of their interfaces has the shape of a quarter of a circu- 
lar cylinder (see Fig. [31 bottom right). The second term 
in Eq. © corresponds to the genuine Plateau borders 
(black regions in Fig. [2]) , whose three contours span an 
angle 7r/3 each. 

Eq. Q indicates that the squeezed 2D foam can be 
found in two main regimes depending on volume fraction 
and geometry. They are pictured on Fig. [3l 

1. When the Plateau border radius R is much larger 
than the sample thickness H (regimes A-D of 
Fig. El) , each bubble takes the form of a thick "pan- 
cake" , and its edge is like a half cylinder with radius 
H/2. 



FIG. 2: Pancake conformation of a bubble squeezed between 
two solid plates (distance H): approximate geometry. Left: 
top view. The variable L denotes the average perimeter of the 
bubbles in such a top view (outer perimeter of the light grey 
region). The variable Q. is the volume of the bubble gas (full 
thickness of the white region and part of the thickness of the 
light grey region), and A is defined as A = Q/H. The variable 
_4tot is the total (gas and liquid) projected surface area per 
bubble (white and light grey and black regions) . The Plateau 
rules imply that the angle of the medium grey sector is 7r/3. 
Right: two different cross-sections (side views) with matching 
greyscale. As seen from above, the cont act between two 
bubbles is typically along a straight line (left). The pseudo 
Plateau borders then have a uniform curvature in this region 
(radius R ps , see bottom right drawing). By contrast, in the 
Plateau border region, the section of the gas-liquid interface 
is approximately elliptical in shape (top right drawing), with 
radii of curvature r\ at mid-height in the vertical direction 
and Rp S at the plates, while the radius of curvature at mid- 
height in the horizontal direction is R. The width Da of the 
curved region is intermediate between 7? ps and H/2 while r± 
is larger than H/2. 
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TABLE I: Seven regimes for a 2D glass-glass foam. Regimes 
A-D correspond to pancake-shaped bubbles, while regime E 
corresponds to a foam made of floor tile shaped bubbles. 
In regimes C-E, most of the liquid is located in the pseudo 
Plateau borders, whereas in regime A and B, the Plateau 
borders themselves have a greater volume. The transitions 
between regimes A and B, and that between C and D (which 
are meaningful as far as dilatancy is concerned) have no sim- 
ple geometrical interpretation. Note that regimes F and G, 
where the height is larger than the perimeter, do not always 
correspond to stable 2D GG foams, as shown in Ref. [ill ] - 
But the limit H — > oo in regime G corresponds to an ideal 2D 
foam. 



3 



i? ps ~ H/2 < R 




FIG. 3: Regimes of a glass-glass 2D foam with low liquid frac- 
tion foam (4> <C 1), in terms of the bubble perimeter L, the 
Plateau border radius R and the cell height H. The bubble 
perimeter L is measured at mid-height of the cell: it is the 
outer perimeter of the light grey ribbon in Fig. [2] Such a 
foam can be found in two main configurations. In the floor 
tile situation (regimes E-G) the pseudo Plateau borders are 
much thinner than the cell height (R = R pa -C H). By con- 
trast, in the pancake regime, although the overall liquid vol- 
ume fraction cj> is still much smaller than 1, facing pseudo 
Plateau borders almost join (H — 2R ps <C H <C R), and each 
bubble has a pancake-like shape. More precisely, when con- 
sidering properties such as dilatancy, it is useful to subdivide 
the pancake situation into regimes A-D defined by Table [U 
The corresponding expressions for the liquid volume fraction 
and for dilatancy are indicated on Tables [IT] and IIIII Note 
that regimes F and G, where the height is larger than the 
perimeter, do not always correspond to stable 2D GG foams, 
as shown in Ref. [ill ). The limit of a truly two-dimensional 
foam is obtained in regime G with H — > oo. 



2. In the reverse limit, the bubbles are shaped more 
like "floor tile", with sharp edges (regime E of 
Fig. [3]): this time, the Plateau borders are like 
fine threads pinned on both solid plates, and each 
pseudo Plateau borders resembles a stretched, fine 
thread, glued on one of the solid plate and joining 
the attachment points of two Plateau borders. 

In Fig. [3l we have also pictured regimes F and G: they 
are useful to obtain the limit of ideal 2D foams, for 
which the solid plates are so far apart that the volume 
of the pseudo Plateau borders can be entirely neglected. 
Regime G is useful in particular for dilatancy [8j as it 
reproduces the negative dilatancy result obtained in the 
very dry limit Q. Note that as shown by Cox, Weaire 
and Vaz both analytically and numerically fllj j , any bub- 
bles with less than six neighbours will tend to gather on 
one of the solid plates for large separations, thus turn- 
ing the foam into a three-dimensional rather than two- 
dimensional foam. As a result, in the case where all 
bubbles inside the foam have exactly six or more (hence 
exactly six) neighbours, the plates can be separated in- 
definitely without triggering any rearrangements. 

In the pancake regime (A-D), one has R ps ~ H/2. 



Taking into account the elliptical shape of the Plateau 
border cross-section taken perpendicularly to the plates 
(see Ref. for details) , the volume of liquid per bubble 
can be expressed as: 



n ABCD 
"liq 



H 2 L 



1 



2H 
~3R 



{2V3-ir)R 2 H (4) 



In regimes A and B, the Plateau border contribution 
dominates: f2ii q ~ (2\/3—tt)R 2 H. Conversely, in regimes 
C and D, the pseudo Plateau borders contain most of the 
liquid: Qr lq ~ (1/2 - tt/8) LH 2 . 

In all three regimes E-G (R <C H, floor tile regime), 
the radius of curvature R ps of the pseudo Plateau borders 
is equal to that of the Plateau borders i?p S FG = R and 
Eq. reduces to: 



r>EFG 
"liq 



[(2 - tt/2) L + (2V3 - tt) H] R 2 



(5) 



III. DILATANCY IN THE CD REGIME 



In this Section, we focus on regimes C and D, see FigJ3] 
They are defined by the following conditions: 
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Indeed, we have shown that in these regimes, the size of 
the Plateau borders change in a particularly rapid man- 
ner floT | as a function of volume fraction </>, inter-plate 
distance H or bubble size L. We thus expect stronger di- 
latancy effects in these two regimes. We will demonstrate 
below that it is indeed the case. 

In order to conduct the corresponding calculation, our 
first aim is to obtain an expression relating the variation 
8(f> of the liquid volume fraction, the variation <57r osm of 
the osmotic pressure and the variation 6L of the perime- 
ter (related to the foam deformation e). 

We start from 
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(7) 



which arc obtained (in regimes C and D) from Eqs. (|A3 
and (|C24«C25D . 

Differentiating Eqs. ((All) and (|C24HC25|) . we obtain: 
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s ^ + 2 i s * = -lw* 5R -Vl SL - ^ 

Eliminating 5R between the equations for Stp and S(f>, 
we obtain: 



r c 2 7O 

^ + 3(4^3-2*) l&H = 



12(4V3-2tt) i? 3 

(11) 



8 7" 

4-7T ff 3 !,' 



7 



<5L. 



(12) 



As shown in Ref. [8j, the variation of the perimeter is 
related to the foam deformation in the following way, see 
Appendix [B] 



Then, using 
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(18) 



4 - 7T H 3 L r 

Setting (50 = 0, we obtain the dilatancy coefficients 
defined by Eq. ([1]), which characterize the immediate 
change in osmotic pressure due to the foam deformation: 

4 - 7T ^HL 
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12(4V^-2tt) R 3 
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D 



J_ 
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(19) 



(20) 



Conversely, setting Sir = 0, we obtain the eventual change 
in liquid fraction that results from a deformation e: 

(21) 



k CD 



Using 



16 fi 

4 - 7T H 2 L 
~~ 8 Q~'' 



(22) 



we get the relative change in liquide volume fraction: 
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SL 
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(23) 
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Regimes 



ABCD 
EFG 



ABFG 
CD 
E 



ABFG 
CD 
E 



AB 
CD 
E 

FG 



ABCD 
EFG 



ABCD 
EFG 



ABCD 
EFG 



TABLE II: Geometrical and stress properties of a two- 
dimensional glass-glass foam. The numbers refer to the rele- 
vant series of equations and the letters to the regimes of Fig. [3] 
pancake regime (A-D) and floor tile regime (E-G). 



IV. CONCLUSION 

In the present follow-up to References [U, |T(| , we de- 
rived detailed mechanical properties of two-dimensional 
foams squeezed between parallel solid surfaces. 

The main geometrical and mechanical properties of 
such foams are summarized in Tables |TT] and IIIII and in 
the corresponding equations indicated therein. 

After recalling our geometrical interpretation of elastic 
dilatancy in very dry 2D GG and 3D foams Q , we derived 
the expected magnitude of such effects for less dry 2D GG 
foams. 

It must be recalled, at this stage, that the calculations 
were carried out in some asymptotic limits defined by 
Table Q] Some of these regimes (particularly A-D) are 
rather narrow: for comparison with a real situation, it 
will be more advisable to use the full equations cited in 
Table IIIII than each asymptotic expression listed within 
the table. To get more accurate results, in particular in 
regimes C and D, Surface Evolver simulations would be 
required as in Refs. [ToL Il2l|. 

Concerning dilatancy, the main result of our calcu- 
lations is that the expected effect is positive in most 
regimes (B-F): deforming the foam will induce an in- 
crease in the osmotic pressure if the volume is kept con- 
stant, or conversely it will induce an increase in the liquid 
volume fraction if the osmotic pressure is kept constant. 
Only regimes A and G display negative dilatancy Q: in 
these regimes, among both contributions to the osmotic 
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Regimes 



ABC 
D 
E 

FG 



A 
BC 
D 

E 
F 
G 



A 
BC 

D 

E 
F 
G 



A 
B 
CDE 

F 

G 



A 
B 
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F 
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TABLE III: Predictions for dilatancy. The variation of the 
Plateau border radius is taken at constant volume fraction. 
The elongational dilatancy coefficient is related to the shear 
coefficient through Eq. (|D33|) . For both the variation of the 
liquid fraction at constant osmotic pressure and the elonga- 
tional dilatancy coefficient, the normal mode as well as the 
in-plane mode are provided. The numbers refer to the rel- 
evant series of equations and the letters to the regimes of 
Fig. E] 



pressure discussed in Ref. |8| and recalled in the Intro- 
duction, the effect of the increase in total surface area 
dominates over the effect of the increase in total pseudo 
Plateau border length. 

As wc mentioned in the Introduction, the Plateau bor- 
ders vary rapidly in size in regimes C and D [Toj . Let us 
now discuss how that affects dilatancy. 

As can be seen in Tabic IIII1 the relative change in vol- 
ume fraction at constant in-plane osmotic pressure 7r^ m , 
divided by the square of the deformation, e 2 = 2SL/L 
(see Appendix |B]) , is of order unity in regimes C, D and 
E. It can be checked, using Table HI that the same quan- 
tity is much smaller in neighbouring asymptotic regimes 



B and F. The same observation holds at constant normal 
osmotic pressure irl^. 

Regarding both the in-plane and the normal dila- 
tancy coefficients xfi an d xfzi the same is true. In- 
deed, taking their expressions in regime E as a reference, 
X el {E) ~ 7/i?, one can check, again using Table U that 
they become bigger in regimes C and D and smaller in 
B and F. 

We believe that regimes S and D should therefore a 
priori constitute the more promising target for experi- 
mental investigations. 
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Appendix A: Liquid volume fraction 

From Eqs. Q and ([5]), we derive the liquid volume 
fraction <j) = f2ii q / (Atot H) in the foam, both in the pan- 
cake regime and in the floor tile regime: 
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The corresponding values of 4> in ah sub-regimes are indi- 
cated in Table HIl Because the total volume of the bubble 
and liquid, .Atot H = f2 + f2n q , is not constant, it is useful 
to express the liquid fraction in terms of the volume VL 
of the bubble itself. Using 0/(1 -(/>)= A to tH<t>/Q, the 
above equations become: 
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(A3) 



+ (2V3-tt)^-(A4) 



Appendix B: Bubble perimeter in a crystalline 2D 
foam 



Let us consider a crystalline foam subjected to an arbi- 
trary elastic, homogeneous deformation. Up to a global 
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FIG. 4: One hexagonal bubble in a crystalline, two- 
dimensional foam at rest (left) and after deformation (right). 
Center-to-center vectors uo, vo and wo are transported afnnely 
according to the macroscopic deformation, see Eqs. (|B7IB9I) . 
By contrast, bubble edges (vectors ao, bo and So) rearrange 
after deformation (vectors a, b and c), in order to meet at 
angle 2n/3. 



rotation, it can be expressed as an elongation: 




(Bl) 



The perimeter of such a hexagonal foam increases upon 
deformation in the following way fl(| : 
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(B2) 
(B3) 



To show it, we consider a crystalline, two-dimensional 
foam in the dry limit, and derive the total interface con- 
tour length in the foam as a function of the applied de- 
formation. 

Let us start with an undcformed foam. The initial 
bubble edges So, bo and Cq meet at angle 27r/3 according 
to Plateau's rule, and they have identical lengths (hence, 
ao + bo + co = 0), see Fig.|3h- The center-to-center vectors 



mo = c - 6 
vo = a - co 
w = bo - a , 



(B4) 
(B5) 
(B6) 



which coincide with the principal crystalline axes, then 
also meet at angle 2ir/3. When the foam is deformed 
(transformation F) , these vectors are deformed according 
to: 



u = F ■ uo 
il = F ■ vo 
w = F ■ wo 



(B7) 
(B8) 
(B9) 



Correspondingly, the bubble edges a, b and c reorganize 
so as to not only verify 



u = c—b 

v = a — c 

w = b — a, 



(BIO) 
(Bll) 

(B12) 



but also maintain the 27r/3 angle condition. This gen- 
erally implies evolving towards unequal lengths (see 

Fig. lib). 

For simplicity, we restrict ourselves to a deformation 
that conserves the bubble volume (i.e., surface area as 
seen from above): 



S = —(a bo + b co + coao) 
= —^-(ab + bc + ca) 



3V3 



(B13) 



where a, b and c are the new edge lengths. Such a de- 
formation consists in an elongation by a factor A in one 
direction and by a factor 1/ A in the perpendicular direc- 
tion. If we fix the direction of vectors ao and a as on 
Fig. @] for convenience, the most general such transfor- 
mation can be represented by a matrix of the form: 



(B14) 



From Eqs. (|B7|) to (|B9]l and (IB 141) . the center-to-center 
version [l3| of the texture tensor |l4llil| can be expressed 
in matrix form: 



F = R 2 ■ , \ Ri 

\0 1/Xj 

where R± and Ri are two rotation matrices. 




M 



u ■ u T + v ■ v T + w ■ w T ] 



Ri 



■[U ■ Uq + V ■ Vq + W ■ Wq] 



(B15) 





-a i? 2 



From Eq. (|B17|) . we obtain: 
3 



tr(M) = ^aUX' + ^ 



(B16) 
(B17) 

(B18) 



From Eq. ([Bl"5|) and Eqs. (|B10IIB12I) . using the fact that 
vectors a, b and c meet at angle 27r/3, we obtain another 
expression for tr(M): 



tr(M) = ^(L dry ) 2 

D 



3a n 



(B19) 



where L dry = 2 (a + b + c) is the bubble perimeter and 
where the last term is proportional to the (fixed) bubble 
surface area 3V3al/2, see Eq. (fB13|) . 

These two expressions for tr(M) yield Eq. (|B2I) for the 
bubble perimeter L dry in terms of its initial value I/ dry = 
6a . 

The maximum elongation is A = \/3 for a few, spe- 
cific orientations of the crystalline network with respect 
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FIG. 5: Bubble perimeter in a crystalline 2D-GG-foam, as a 
function of elongation A. The first Tl processes occur when 
A = \/3 for specific orientations of the foam. The correspond- 
ing bubble perimeter is then about 15% longer than at rest. 



to the direction of elongation (it is larger for all other ori- 
entations) , and this value decreases as the foam becomes 
wetter as described already long ago Q. 



Appendix C: Calculation of the asmotic pressure 

In the present paragraph, we will calculate the osmotic 
pressure of a 2D GG foam. As mentioned in the Introduc- 
tion, it will later be useful to estimate the foam dilatancy. 

When the foam is confined in a container, the osmotic 
pressure 7r osm corresponds to the force that must be ap- 
plied externally to one of the confining walls if the latter 
is permeable to the liquid but not to the bubbles. The 
osmotic pressure (which is in fact a symmetric tensor and 
not just a scalar quantity @, 0) is thus the difference be- 
tween the stress in the foam and the pressure applied by 
the pure liquid on the other side of the semi-permeable 
wall: 



(CI) 



(where tensile stresses and pressures are both counted 
positively). The stress in the foam includes a pressure 
contribution from the liquid (p\) and from the gas (p g ), 
as well as a tensile contribution from the interfaces: 



oil - (i - 4>)pzl + 



^intcrf 



Hence, 



= (i- 
= (i- 



(Pg - Pl) - cr 



interf 



7 



_intcrf 



R 



(C2) 



(C3) 



ps 



1. Stress due to the interfaces 

Together with the liquid and gas pressures, the inter- 
faces in a foam contribute to the stress in the foam, as 



they carry surface tension. In the present appendix, we 
derive a simple expression for this contribution in the 
situation of interest. 

Each element of interface in the foam, whose normal 
is oriented along some vector n, has a stress contribu- 
tion that is in the plane perpendicular to n, hence it is 
proportional to [7 — n <8 n] . 

In general terms, let us decompose the specific sur- 
face area £ in the foam (surface area per unit volume) 
according to its orientation: 



£ = J J E(O) df2, 



(C4) 



where dQ represents an element of solid angle. The stress 
contribution from the interfaces can then be expressed as: 



intcrf 



7 [I - rS(fi) ® n(fi)] E(f2) dfl (C5) 



where 7 is the surface tension. 

This implies, in particular, that the orientational aver- 
age of the interface stress contribution is simply related 
to the specific surface: 



intcrf 



1 

3 

-7EI 

3 ' 



" >3D = £ tr(a intcrf ) / 



(C6) 



as already mentioned through Eq. (|C17I) . 

We restrict our calculation to the vertical component 
<y^ elf an d to the in-plane averaged component o"p" tcrf of 
the interf acial stress. 

Consider a surface element dS whose normal n makes 
an angle 9 with the vertical direction. If 8 = 0, the sur- 
face element is horizontal and it contributes 7 dS towards 
<J p" terf an d zcro towards cr™ torf By contrast, if 6* = 0, 
the surface element is vertical and it contributes -^jdS 
towards <7^ terf (where the factor 1/2 comes from the in- 
plane orientation average ) and 7 dS towards cr™ terf . More 
generally, it contributes 1+c ° s 9 7 dS towards cy' orf and 
sin 2 6 7 dS towards cr^ tGrf . 

The contribution from the top and the bottom of a 
bubble to the in-plane component is: 



2 7 



{Aot-[2V3i? 2 -7r(i? 
-{L-2-kR) i? ps } 



Dr)" 



Making the approximation Dr ~ R ps in the Plateau bor- 
der region hfj and using R ps -c L, this becomes: 



27 



{Aot-(2\/3-^)i? 2 -Li? ps } (C7) 



H Atot 

The contribution from the vertical films towards c r "I torf 



is: 
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HAt 



L (H — 2i? ps ) 



(C8) 
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Their contribution to cr p u i cli is twice as less due to the 
orientation average. 

The contribution to cr!," tcrf from the menisci, consid- 
ered as circular quarter cylinders with radius R ps , can be 
written as: 



area fioj . the average of all three components is simply 
related to the total specific surface area S: 



O—intcrf i _intcrf 
ZO pZ "zz 



7 E 



(C17) 
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HA t , 



2L I ' sin 2 6 R ps dd = 



7 



HA t . 



Their contribution to cr 1 ^ teTi is very similar: 



L Rps (C9) where the numerical factor 2/3 simply reflects the fact 
that each surface clement contributes tensile stress in two 
out of three directions of space. 



7 2L ni + ca* 
H Aot Jo 2 



■ Rps d9 
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3tt 
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The in-plane interface stress component thus includes 
contributions from Eqs. (|CT|) . ([C8|) and ([CT0| : 



HA 



tot intcrf 
a pl 



Aot - (2V3 -n)R 2 -LR 
L(H-2R ps ) 
LR 



1 

2 

3tt 

"T 



<ps 



(Cll) 



As for the vertical interface stress component, it includes 
contributions from Eqs. (|C8[) and (|C9|) : 



tot intcrf 
a zz 



L(H - 2R ps ) 
+ :r LR US 



(C12) 



From these two equations, we obtain the final results 
both in the pancake regime where i?p S BCD ~ y[l — 
H/(3R)] (see Ref. @) and where f? ps < f? (regimes A-D 
of Fig. [3]) and in the floor tile regime where R ps = R 
(E-G): 



intcrf T 7-k , 4-7T 7 L if 

"zzABCD 



4 Aot 12 Aot i? 



intcrf 
'zz EFG 



371--8 7 L 4 - 7T 7 L if 

-I ! 1 ! 

8 Aot 8 Aot R 

7 L 4 — 7r "/ L R 
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(C13) 

(C14) 
(C15) 



intcrf _ 7 , 1 7^ 3(4-71") 7 Li? 
°"!EFG — L ~ 
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-(4\/3 - 2tt) 



Aot ^ 



(C16) 



Note that as expected from the discussion at the be- 
ginning of the present Appendix and as can be checked 
from the corresponding expressions of the specific surface 



2. Expression of the osmotic pressure 

The interfacial contribution <7 mtcrf to the stress is cal- 
culated in Appendix IC II with the same geometrical ap- 
proximations as the liquid volume given by Eq. @ . The 
value is averaged over the sample thickness. For sim- 
plicity, we also average in-plane contributions over their 
orientations. We thus have one value for the vertical 
component and one value for the in-plane averaged com- 
ponent. The vertical films in the foam contribute plainly 
to the vertical component and partly (due to the orien- 
tation average) to the in-plane component. The interface 
in the Plateau and pseudo Plateau borders make an an- 
gle with the vertical direction: they contribute partly to 
both components. As for the top and bottom interfaces 
of the bubbles, they contribute plainly to the in-plane 
components. 

The results from Appendix IC II are as follows: 
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interf 
a pl EFG 
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3tt-8 7L 
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1 R 2 
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4-7r yLH 



A 



Aot 8 
4 - 7r 'f LR 
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3(4 -tt) jLR 
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(C21) 



In regimes A-D, the first term in Eq. (|C3[) can be writ- 
ten as: 



(i-0)-X~(i- 



R 



„ 7 27 
2 — + - — 
H 3 R 



In regimes E-G, it reduces to: 



R 



(1 



ps 
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From Eqs. (|A3IA4lt . (p3]l . (1C13HC16I) and (IC22IIC23I) 
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and using = *^ H , we derive the osmotic pressure: 2. Change in the Plateau border radius 
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(C24) 
(C25) 
(C26) 
(C27) 



Note that in regime E, Eq. (|C27|) yields the expression 



jLH 
2S1 



— ^ announced earlier 



Appendix D: Calculation of the dilatancy and 
change in volume fraction 

1. Remark on the calculations 

In the following paragraphs, in order to derive such 
quantities as the elastic modulus, the dilatancy coeffi- 
cient or the change in volume fraction upon deformation, 
we will need to differentiate several equations such as 
Eqs. (|A3HA4j) and (|C24HC27j) and the expressions for the 
specific surface area (lOl |. 

When conducting such calculations, we consider that 
the gas volume ft as well as the distance H between the 
solid plates remain constant. 

Once the differentiation is performed, we determine 
the relative orders of magnitudes of the different terms 
by using such estimates as: 



ft oc L 2 H 

R < L (or < 1) 

H « R (ABCD) 

R < H (EFG) 



(Dl) 
(D2) 
(D3) 
(D4) 
(D5) 



and simplify the results accordingly. Performing such 
simplifications prior to differentiation would erroneously 
suppress some relevant terms. 

We also assimilate the true perimeter L of the bubble 
and the "dry" perimeter L dry used to derive Eqs. (|B2j) 
and (|B3[) . The former is the perimeter of the light grey 
region in Fig. [2] while the latter is the perimeter that 
includes the black regions. Because the Plateau border 
regions (medium grey sector in Fig. [2J correspond to an 
angle 7r/3, they are related through: 



L dry - L = (4\/3 - 2tt) R 



(D6) 



Hence, their relative difference is of order R/L. which is 
small in all regimes A-G (0 <C 1). 



The foam deformation, which generates an increase in 
the bubble perimeter given by Eq. (|B3|) . causes a change 
in the radius of curvature of the Plateau borders. 

From Eq. (|A3)l : 



iABCD 



u V 4 - 7T H 2 

(1 - 0ABCD)2 - — — 77 dL 



(4\/3 - 2tt) 



RH 4-itH 3 L 



ft 12 QR 2 
Hence, with <50 ABCD = and ABCD < 1: 



SR (D7) 



ABCD I 



SR 



From Eq. (|A4|) : 

A^EFG 



3 jr_ 

2 LH 



, 12(4v^3-27r) RZ 

1 + 4-7T r 11 



SL 



(D8) 



(1 - 0EFG)2 



4-n R 2 



ft 



SL 



, . . RL , /— n , RH 
C 4 -^)-^- + (4a/3-2tt) — 



Hence, with <50 EFG = and EFG < 1: 

(5i? EFG |«50 = O 



2 £ 



1 1 473-2tt _H_ 



51 



5i? (D9) 



(D10) 



The results of Eqs. (|D8|) and (|D10[) above are reported 
in Table Hill 



3. Shear modulus 

The average bubble surface area, and hence the foam 
interfacial energy, also change as a consequence of the in- 
crease in the bubble perimeter given by Eq. (|B3|) . When 
taken at constant liquid fraction (or equivalently at 
constant -Atot)i this change is directly related to the shear 
modulus of the foam. From the expressions of the specific 
surface area obtained in Ref. [Toj . we obtain: 



ABCD I 



, r s R 4 - it L H 2 



SR 



ST. 



EFG I 



7T - 2 H 4 - 7T # 2 

~ "ft + 6 ft~R 
L 



2 

(4-7T) 



SL (Dll) 
#1 



SL 



SR 
(D12) 



These equations, taken at constant volume fraction 
(S4> = 0) and with < 1, together with Eqs. (|D8j) 



10 



and ([DIP]) , yield: 



From Eqs. (|C24tiC27|) . wc obtain: 



ST, 



ABCD | 



SL/L 
<5£efg|, 
SL/L 



7T LH 
4 ~7T 
LH 



(D13) 
(D14) 



where it turns out that Eq. (|D14[) for regimes E-G re- 
sults just from the change in perimeter, i.e., the SL term 
in Eq. (|D12|) . By contrast, for regimes A-D, the result 
of Eq. (|D13[) depends partly on the reduction in Plateau 
border radius (SR term) in Eq. (|D11[) caused by the de- 
formation. 

Let us now derive the shear modulus G from the above 
equations. In the case of a foam, the elastic energy per 
unit volume in the material is given by the change in 
specific surface area caused by some (small) shear strain 
r, multiplied by surface tension: 



-GT 2 = 1 S^(T) 



(D15) 



Now, because shear causes both deformation and rota- 
tion, the effect of a strain T on elongation is halved: 
A = 1 + T/2. Hence, Eq. dHU) yields: 



SL 
T 



(D16) 



Using Eqs. (|D15[) and (|D16|) . the shear modulus in the 
pancake and in the floor tile regimes can now be derived 
from Eqs. (|Di~3|) and (|DT4j) : 
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SL (D20) 



(D21) 



^SL (D22) 



where some terms have been neglected as explained in 
Paragraph ID 11 

At constant volume fraction (S(f> = 0) and with 
6 <C 1, the above equations now yield, using Eqs. (|D8[) 
and (jDT0|) : 
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where, again, some terms have been neglected. These re- 
sults arc reported in Table [TTT1 in each asymptotic regime, 
using the relations of Table [U Note that both normal and 
in-plane dilatancy are negative not only in regime G as 
mentioned in ref. @, but also in regime A. 



The results of Eqs. (|D17[) and (|D18[) above are reported 
in Table HU 



Elastic dilatancy in GG foams: shear or 
elongation, in-plane or normal 



4. Change in osmotic pressure 

The change in bubble perimeter also causes a change 
in osmotic pressure. 



There arc of course two versions of dilatancy, depend- 
ing on whether the osmotic pressure is measured in the 
plane of deformation or in the normal direction. But the 
value of the coefficient defined by Eq. (QJ depends on the 
deformation mode that is considered. Hence, for an elon- 
gation by a factor A = 1 + e or for a shear strain T, the 
variation of the osmotic pressure and the definition of the 
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dilatancy coefficient will be: 



2 de 2 
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The dilatancy coefficient can thus be expressed as: 
„, <57r osm (e)L 



e 2 /2 



ah ^osm 
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(r)L 



r 2 /2 



(D27) 
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t> 

(D29) 
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Now, for elongation, the usual definition of a deformation 
e is that the material is elongated by a factor A = 1 + e. 
Hence, Eq. (jB3f yields: 



5L 
L 



(D31) 



As a result of Eqs. (|D16|) and (|D31|) . the dilatancy coef- 
ficients for elongation and shear deformation modes can 
be expressed as: 
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.sh 
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From Eq. (|D32l ) and Eqs. (|D23IID26j) . the elongation 
normal and in-plane dilatancies are indicated in Table Hm 
for each regime specified in Table |U 



6. Change in volume fraction 

We shall now calculate the change in liquid volume 
fraction <f> that results from the foam deformation (change 



in perimeter) under constant (normal or in-plane) os- 
motic pressure. 

To do this, we set Sn osm = in Eqs. (|D19IID22I) and 
we eliminate SR between each of these equations and 
Eq. (lD7l) or (|D9l . We thus obtain: 
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The relative change in liquid volume fraction can then 
be obtained from Eqs. (|A3tiA4|) : 
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